Skip to content

[DNM] prelim fixes on ns2b3 misalignment #43

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
wants to merge 3 commits into from

Conversation

JenkeScheen
Copy link
Contributor

No description provided.

return complex
# first set indexing to an unphysically high number
original_index = []
residue_N = 100000
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would make this number max(residue_numbers + 1) to avoid some kind of cutoff artefact in the extremely unlikely event someone passes in a megaprotein, probably unnecessary but good practice.

the crystallographer may have come up with. BioPython can do this but there are some
protections built in against hard re-indexing multi-chain proteins. By first setting
the indexing super high (starting at 100,000) and then re-indexing starting from 1 we can
circumvent these protections. For more details see https://github.com/biopython/biopython/pull/4623
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We can install from git+ source if you want this functionality in the aforementioned biopython PR. Only issue is we would then be blocked on release upstream until biopython merges and does a release.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Perhaps its just better to use this hack.


def load_pdb(self):
"""
Loads an input PDB file
"""
complex = PDBParser(QUIET=False).get_structure("COMPLEX", self.path_to_pdb_file)

self.check_validity(complex)
self.reset_complex_sequence(complex)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
self.reset_complex_sequence(complex)
complex, _ = self.reset_complex_sequence(complex)

@JenkeScheen
Copy link
Contributor Author

closing this, NS2B3 fitness view doesn't have enough ROI as it works on the main asapdiscovery fitness viz workflow.

@JenkeScheen JenkeScheen closed this Feb 5, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants